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INTRODUCTION 

The  following  is  a  description  of  a  program  for  a  hand-held  programnaPle 
calculator  which,  in  effect,  interpolates  values  of  temperature  from  the 
calibration  table  of  a  silicon  diode  thermometer  over  a  full  4  to  300K 
temperature  range.  When  making  temperature  measurements  in  the  laboratory, 
the  voltace  across  a  forward  biased  calibrated  diode  is  read  from  a  digital 
voltmeter.  With  this  program,  the  voltage  reading  is  keyed  into  the  calculator 
and  the  interpolated  temoerature  is  calculated  quickly  and  accurately.  The 
computational  technique  oresented  here  is  superior  to  techniques  generally 
used  since  the  data  reduction  scheme  allows  more  information  to  be  stored, 
allows  for  a  simple  calculator  routine  resulting  in  very  short  computation 
times  and  results  in  an  interpolation  having  the  accuracy  of  the  original 

data  table.  The  scheme  is  perfectly  general  and  cjo  oe  used  to  interpolate 

0 

any  data  table. 
ucNEML  SCHEME 

A  camion  recurring  proolem  in  laboratory  work  is  the  necessity  of 
interpolation  f'-om  a  table  of  values.  An  example  is  interpolation  **-om  a 
calibration  table  furnished  with  a  cryogenic  silicon  diode  thermometer  where 
pairs  of  temperature  and  voltage  values  are  given/1  During  an  experiment, 
the  diode  voltage  is  read  from  a  digital  voltmeter,  then  the  corresponding 
temce^ature  value  must  be  interpolated  from  the  calibration  table.  It  would 
be  very  convenient  to  be  able  to  make  this  interpolation  quickly  and  with  a  high 
degree  of  accuracy  while  the  experiment  Is  in  progress.  A  hand-held  programmable 
calculator.  In  our  case  a  Texas  Instruments  TI-59,  can  be  used  to  accomplish 
this. 
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Before  this  is  done,  a  number  of  problems  must  be  overcome.  The 
calibration  table  Supplied  with  a  diode  typically  may  have  90  or  more  pairs 
of  temperature  and  voltage  values.  This  would  exceed  the  maximum  available 
memory  capacity  of  the  calculator  which  is  about  100  locations.  Interpolation 
can  be  accomplished  by  the  evaluation  of  a  least  square  fit  polynomial, 
nowever  the  calculator  is  too  slow  to  do  other  than  the  simplest  fitting 
routines  in  a  short  period  of  time.  An  attempt  to  fit  a  polynomial  equation 
to  the  entire  A  to  300K  calibration  using  a  large  computer  did  not  result 
in  the  desired  accuracy,  "he  following  procedure  was  subsequently  devised. 

"he  values  from  the  calibration  table  were  divided  into  overlappinn 
temperature  intervals.  Then  for  each  Interval,  a  least  square  polynomial 

i'?' 

*1t  was  made  using  a  Forsythe  fitting  method'"'  programmed  on  an  IBM  370 
compute'*.  Any  least  square  polynomial  *1 1  with  sufficient  accuracy  is  adequate. 
~he  ’east  square  polynomial  #it  program  gives  an  ms  error  for  each  polynomial. 
"He  temperature  intervals  were  adjusted  by  trial  and  error  in  order  to  obtain 
a  satisfactory  rms  error  over  the  *ul  1  A  to  30CK  span.  In  the  case  considered 
nere.  it  was  *ound  that  1C  polynomials  each  of  Ath  order  were  needed  and  that 
tKe  length  pf  the  temperature  intervals  varied  over  the  calibration  range. 

All  intervals  have  a  one  point  overlap  with  the  adjoining  intervals.  Using 
the  coefficients  found  to  16  place  precision  by  the  IBM  370  computer,  the  rms 
er»-or,  dg^Q,  of  each  polynomial  fit  is  given  in  Table  I.  The  accuracy  obtained 
preserves  the  accuracy  of  the  calibration  table,  "he  full  calibration  table 
was  thus  reduced  to  12  sets  of  polynomial  coefficients,  a0>  a1 ,  a^.  a^,  a4, 
which  could  all  be  stored  In  the  calculator.  A  relatively  simple  program  could 
then  be  used  to  calculate  any  temperature.  T,  between  4  and  300K  by  evaluating 
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T  *  a^  V  ♦  a ^  V  «■  a?  V  ♦  V  ♦  aQ,  where  V  is  the  voltage  reading  across 
the  silicon  diode  thermometer.  Both  the  coefficients  and  the  calculator 
program  are  stored  on  magnetic  cards  and  a  set  of  magnetic  cards  can  be  used 
to  store  the  coefficients  for  several  different  calibrated  temperature 
sensors . 

Table  I.  The  temoerature  intervals  and  the  corresponding  rms  errors  of  the 
4th  order  polynomial  fit  using  coefficients  having  16  digit  precision  on  an 
IBM  370,  rjjg,  and  coefficients  having  10  digit  precision  on  a  TI-59,  059. 


Polynomial 

• 

Temp.  Interval 

00 

'370 

(«•) 

J59 

(K«) 

1 

4  -  6.5 

1 x 1 0” * 

1 x 1 0~ 1 

2 

6.5  -  9.5 

2x 10“* 

2x10"* 

3 

9.5  -  13 

9x10"' 2 

3x10" 5 

4 

13  -  17 

6x1  O’1 2 

IxlO'5 

5 

17  -  21 

4xl0',: 

2x10’* 

6 

21  -  25 

6xl0"12 

1x10“ 

7 

25  -  32 

7x10'* 

7xlO"5 

3 

32  -  55 

lxl 0~ 2 

IxlO"1 

9 

55  -  35 

IxlO"1 

IxlO"2 

10 

o5  -  155 

3x10"’ 

9x10" 5 

11 

155  -  245 

9x10*’ 

IxlO"1 

12 

245  -  300 

2x10"* 

2x10"’ 
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CALCULATOR  PROGRAM 


The  12  polynomials  used  each  have  five  coefficients  therefore,  60  memory 
locations  are  needed  for  coefficient  storaqe.  The  program  itself  uses  13 
memory  locations  for  the  voltages  which  define  the  temperature  intervals, 

4  locations  for  accimulation  registers  and  one  location  for  the  input  voltage. 

The  program  itself  requires  106  program  steps.  There  are  13  memory  locations 
and  133  program  steps  not  used.  The  number  of  memory  locations  and  program 
steps  can  be  reallocated  so  that  a  maximum  of  fifteen  4th  order  polynomials 
can  be  used  if  necessary. 

The  program  action  is  as  follows: 

a)  The  voltage  drop  across  the  silicon  diode  is  keyed  in. 

b)  The  program  finds  the  corresponding  voltage  interval. 

c)  The  polynomial  coefficients  corresponding  to  this  interval 
are  located. 

d)  The  fourth  order  oolynomial  is  evaluated  using  these  coefficients 
and  the  resulting  interpolated  temperature  is  displayed  rounded 
to  three  decimal  places. 

e)  If  a  voltage  is  entered  which  lies  outside  of  the  range  of  the 
calibration  the  display  will  flash  with  the  value  of  the  voltage 
which  defines  that  calibration  limit. 

Using  the  TI-59  with  the  coefficients  stored  to  the  10  digit  precision  of 
the  calculator,  the  temperatures  were  calculated  at  the  known  voltage  calibration 
points  and  checked  against  the  given  temperatures  from  the  calibration  table. 

This  ms  error,  is  given  in  Table  I  for  comparison.  The  calculation  still 
preserves  the  accuracy  of  the  calibration  table. 

The  time  required  for  the  calculator  to  compute  the  temperature  is  ten 
seconds  or  less.  Thus  it  is  possible  to  take  voltage  readings  and  calculate 
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the  temperature  while  an  experiment  is  in  progress.  A  copy  of  the  TI  program 
is  given  in  the  Appendix  as  Program  I.  Familiarity  with  programming  language 
and  procedures  on  the  TI - 59  are  assumed. 

EXAMPLE 

First,  the  program  must  be  entered  into  the  calculator.  The  intial  step 
is  to  partition  the  calculator's  memory.  This  is  accomplished  by  keying  In 
9  2nd  OP  17.  The  calculator  responds  by  displaying  239.39  indicating  the 
correct  partition.  The  display  is  cleared  to  0  and  the  four  sides  of  the  two 
cards  are  ’*ead  into  the  calculator  (''he  display  must  be  cleared  to  zero  before 
each  card  is  entered).  The  calculator  is  now  ready. 

Assume  a  voltage  reading  of  2.0153  volts.  This  value  is  keyed  into  the 
calculator  and  the  program  is  run  by  pressino  ’A",  "he  calculator  then  ^plays 
tie  corr«spondi ng  temoerature  9.753  npunded  to  3  decimal  places.  The  accuracy 
of  the  result  depends  on  the  accuracy  of  the  original  calibration  and  the 
accuracy  of  the  polynomial  fit  as  discussed  above.  Table  13  in  the  Appendix 
shows  the  values  stored  in  the  memory  locations  after  this  example  was  run. 
''emory  location  00  stores  the  input  voltage  reading.  Locations  1-13  contain 
the  voltages  which  define  the  temperature  intervals  starting  with  the  lowest 
temperature  Interval  first.  Location  14  accumulates  a  value  which  determines 
which  set  of  coefficients  shall  be  used.  Location  20  through  79  contains  the 
"“officients  of  the  polynomials  starting  with  the  lowest  temperature  interval 
first  and  stored  in  the  order;  afl,  a^ ,  a^,  a3>  a^.  Location  82  is  a  counter 
which  equals  the  memory  location  of  each  coefficient  to  be  used.  Location  93 
stores  the  result  of  the  evaluated  polynomial  and  location  34  contains  this 
result  rounded  to  three  decimal  places  which  is  the  value  to  be  displayed. 
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ANOTHER  program 


The  storage  capacity  limits  one  to  using  a  maximum  of  fifteen  4th  order 
polynomials.  The  data  set  may  be  such  that  the  desired  accuracy  cannot  be 
obtained  without  exceeding  the  storage  capacity  of  the  calculator.  This 
problem  occurred  with  one  of  our  thermometers  and  the  solution  used  is 
presented  here. 

It  was  found  that  in  three  temperature  intervals  the  fourth  order  polynomial 
did  not  give  the  desired  accuracy.  So  a  fifth  order  polynomial  was  calculated 
'or  these  cases.  A  program  was  written  for  the  calculator  that  selects  the 
apcroonate  numoer  of  coefficients  and  makes  the  decision  whether  to  use  the 
4th  or  3th  order  equation.  A'ter  selection,  the  polynomial  is  evaluated  and 
the  temperature  displayed,  "he  program,  the  contents  of  the  memory  locations 
and  a  'low  chart  i-m  given  in  the  Appendix  under  Program  II. 

It  should  be  noted  that  this  procedure  (chanqing  some  polynomials  to 
5th  order'  was  an  original  attack  on  the  problem.  Our  present  philosophy 
o'  attack  would  be  to  adjust  the  range  of  the  polynomials  by  trial  and  error 
with  the  least  square  fit  program  on  the  IBM  370  to  give  fits  to  accepted 
limits  jsinq  polynomials  all  of  the  sane  order.  Only  if  this  procedure  fails 
should  a  mixture  of  orders  be  used. 

SOWARY 

A  method  has  been  found  to  interpolate  values  from  a  large  table  of  data 
using  a  hand-held  programmable  calculator.  The  method  uses  a  series  of 
overlapping  least  squared  polynomial  fits,  calculated  by  a  large  computer, 
to  reduce  the  information  contained  in  the  table  to  a  small  set  of  polynomial 
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coefficients.  The  entire  set  of  coefficients  can  be  stored  in  the  calculator. 
Interpolation  is  accomplished  quickly  and  accurately  by  the  programmed 
evaluation  of  the  polynomial.  The  method  is  perfectly  general  and,  in 
particular,  is  ideally  suited  to  the  interpolation  of  calibration  tables 
in  the  laboratory  while  an  experiment  is  in  progress. 


REFERENCES 

1.  Also,  the  raw  calibration  data  is  often  available  from  the 
supplier  upon  request  if  a  temperature  calibration  table  is 
ordered. 

2.  Louis  G.  Kelly,  "Handbook  of  Numerical  Methods  and  Applications', 
Addison-wesley  Publishing  Co.  1967,  pg.  63. 
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USER  INSTRUCTIONS 


STEP 

PROCEDURE 

ENTER 

PRESS 

DISPLAY 

!  i 

Partitioning 

9 

2nd 

OP 

17 

239.89 

7 

c 

Enter  program  and  coefficients  from  cards 
Enter  sides  1,  2,  3,  and  4  (not  necessary 
to  be  in  that  order)  clear  display  to  0 
before  each  entry. 

CLR 

o 

-  —  ] 

1  3 

Enter  Voltage  (V) 

V 

V 

1  4 

Calculate  ’emperature  (T) 

A 

j* 

•If  a  voltage  is  entered  which  lies  outside 

of  the  range  of  the  calibration,  the 

display  will  flash  with  the  value  of  the 

i 

voltage  which  defines  the  calibration 

limit. 

] 

r 

[ 

. 

| 

—  _  .  ,  _ _ _ , 

_ 

USER  OEFINfcO  *EYS  DATA  REGISTERS  (  •* 


LABELS (Op  08) 


W3_S_3C_.SS  _S)_:Z 
ir  _ai  _  as  _  m -.«■ 

(L _  '  _ .  )  _  ^  —  10 S  — .  X  ] 

m  _:=:_se  _x_sa.2: 


Program  I 

TITLE  Silicon  Oiode  Thermometer *9343  PAGE __1 _ 0F  2  Tl  PrOQrQnnmQble 

PROGRAMMER  Carl  G.  Gardner _ DATE  ^3  April  79  PfOQrQm  ReCOfd 

Partitioning  (Op  17)  1 2 ,3  9  8  9|  Library  Module _ Pnnter _ Cards _ 


Cards. 


PROGRAM  DESCRIPTION 


This  silicon  diode  thermometer  interpolation  program  is  designed  to  calculate  the 
temperature  (in  deg.  Kelvin)  correspondi ng  to  a  given  voltage.  The  calculator  performs 
a  search  routine  to  find  the  correct  coefficients,  then  calculates  T  using  the 
predetermined  coefficients  of  a  least  square  polynomial  fit. 

T  -  a4  V4  ♦  a3  V3  ♦  a2  V2  ♦  a,  V1  *  aQ 

Range  of  validity  4.0  to  300K. 
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KEY  COMMENTS 


COMMENTS 


CHECK 


RECALL 


N7ERVAL 


RECAL 


CLEAR  7 


ROUND  OFF 


IN7ERP0LA7ED 
'E-MPERATURE  ~ 


MERGED  COOES 

ugg  a  n  sh-  o 

uq  a  -3  «n  r*i  Mnn 

MOO  '*  v<i  a  97  -n  ^ 


Program  I 

tttlf  Silicon  Oiode  Thermotneter  *9343 
PROGRAMMER  Carl  G.  Gardner _ 


of— 1—  71  Rogrommobb 
date  23  APril  79  Coding  Form 


LOC  c 
~000 
'  '001 
’]002 
H003 
.  4 004 
i  005 
(006 
'’007 
008 
00'? 
70 1  o 
Ion 
;  012 
[  '0  1  3 
C  14 


U-5 
-038 
L  0  37 
[-038 
C '  ■- 
04  0 
I  '-4  1 
_  042 
;  ^ : 

,  044 
045 


COMMENTS  ON  PROGRAM  I 


STEPS 
00  to  13 

14  to  28 

29  -  38 

39  -  91 


92  to  106 


Pressing  Label  A  starts  program,  stores  voltage,  V]  for 
wnich  temperature  is  to  be  calculated  into  Memory  Location 
00  and  makes  comparison  to  insure  V  is  within  correct 
range  (in  this  case  2.3822  >  V  >  .33345). 


Locates  correct  voltage  range  by  comparing  V  to  hign  value 
of  voltage  in  each  of  12  ranges.  (These  values  are  stored 
in  Memory  Locations  01  to  13).  Locator  Register  (Memory 
Location  14)  is  initially  set  equal  to  1. 


Calculates  the  starting  Memory  Location  of  the  coefficients 
of  the  polynomial  selected.  Stores  starting  location  in 
"emory  Location  32. 


Evaluates  the  4th  order  polynomial  equation.  The  value 
correspondi ng  to  the  location  of  the  first  polynomial 
coefficient,  a0,  is  in  Memory  Location  32.  'he  next 
coefficient,  aj,  is  found  by  Incrementing  Memory  Location 
32  by  1.  As  the  values  of  the  terns  of  the  Dolynomial 
are  calculated,  they  are  summed  into  Memory  Location  83: 
When  finished  evaluating  Memory  Location  83  holds  the 
answer. 


Rounds  the  value  in  Memory  Location  S3  to  3  decimal  places, 
stores  the  rounded  value  in  Memory  Location  34,  clears  the 
fixing  condition  back  to  full  size  and  displays  the  rounded 
value  from  ‘lemony  Location  34. 
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Table  II  *  Memory  Locations  for  Program  I  after  an  Interpolation  calculation 
r‘or  2.0153  volts. 


MEMORY 

location 


00 


02 

03 

04 

H5 


VALLE 


0  1  f  2 


202  1 
0267 
SO?' 

L.L.~C' 


MEMORY 

LOCATION 


value 


45 

46 
4' 
4  S 
40 


2074. 966311 
-6208. 5364x4 
7116. 287458 
-3654. 244386 
707. 3400023 
404993. 1391 
-  i  396680.  6 32 
1806751.  5  9  9 

?95x.  ‘■‘"i 


_  i  -  :  :  5c 

224076 


c  - 


609130. : 55 
-230  3  606.  3  S-7 
3232707. 3  ~ 
-20750 35.  3 4~ 
491881. 7346 

-  4  C  1  1  .  4  ;  J  - 

1  92639.  31-2 

-^?~*3*2.r^  c  '  _  3. 

- 4 ' 5 9  42-2^ 

254. 7583636 


Program  II 

nn.Ei.UJwn  Jigflg  Theraoffie;er _ 

Interpolation 

PROGRAMMER  Ellen  Fyffe _ 

Partitioning  (Op  i 7)  J  .9,8  .9  ;  Library  Module  . 


.  PAGE  _J _ OF  3 

DATE  1 0-6-73 


Tl  Progrommoble 

Program  Record 

_ Pnnter _ Caros  _ 


PROGRAM  DESCRIPTION 


This  silicon  diode  temperature  sensor  interpolation  program  is  designed  to 
calculate  the  temperature  (in  degrees  Kelvin)  corresponding  to  a  given  voltage.  The 
calculator  performs  a  search  routine  to  find  the  correct  coefficients ,  then  calculates 
J  using  the  predetermined  coefficients  of  a  least  square  polynomial  fit. 

T  •  a4  V4  ♦  33  V3  «-  32  V2  ♦  aj  V  ♦  ag 


(NOTE :  ‘lot  valid  for  voltages  greater  than  2.33785  volts,  (1.9K).  Range  of  validity 


1  Partitioning 

2  Enter  voltage  (V) 

3  Calculate  temperature  (T) 


Program  II 

title  Silicon  Diode  Thermometer 
Interpolation 

PROGRAMMER  Ellen  M.  Fyffe 


LOC  COOG  KE> 

bno'  'miBtr 
om--n  « 


002 

42 

STO- 

t)03 

00 

00 

1004 

32 

x:t 

^05 

01 

1 

A 

006 

42 

STD 

*007 

14 

14 

008 

73 

RC  ♦ 

009 

14 

14 

010 

3  *> 

I  MV 

Oil 

I  ( 

GE 

012 

00 

00 

013 

20 

20 

014 

01 

4 

013 

44 

: ;  M 

016 

1  4 

.  4 

01' 

6  1 

j  » 

0 1  i 

oo 

■ 

019 

08 

os 

020 

0" 

021 

^  A 

•  ^ 

010 

4  3 

r r  l 

go: 

1  4 

«• 

;C14 

*• 

j  I. 

,023 

00 

00 

02* 

02” 

6 1 

3"2 

028 

00 

00 

029 

•4 

4  ] 

,030 

•• 

- :  _ 

1031 

1  •* 

1 4 

032 

„  c 

033 

nc 

c 

Q2  a 

©  c 

+ 

033 

01 

4 

* 

0  36 

02 

0 

03' 

35 

s 

!038 

4  _ 

:  T2 

0  *  Ci 

;040 

61 

GTQ 

04 1 

00 

00 

042 

•4  2 

;  4 : 

03 

tr 

044 

^  •> 

2>  *1 

’  *  T 

043 

4  3 

z 

•  •— 

046 

1  • 

* 

4  * 

L  + 

04  ~ 

J  — 

*1 4 

00 

00 

04  9 

e  t* 

e  “ 

0  30 

£  1 

G'D 

f  *!  4 

a 

•  j* 

:c: 

■* 

■  c  : 

4  ‘ 

s  *  _ 

COMMENTS  I  LOC 

_  035/^65 

_ 03£-~05 


2  OF  3 
10-6-78 


I  COMMENTS 


T1  ftogrommoble 

Coding  Form 


zt=z= 


1  • 
1 


060 

93 

s 

061 

42 

STO 

062 

82 

82 

063 

29 

CP 

064 

06 

6 

06^ 

A  A 
0.1 

x:t 

06 

43 

RCL 

0*.’ 

14 

14 

068 

6  • 

EQ 

069 

01 

01 

070 

49 

4  9 

0^1 

0  —  *> 

i 

L»  1  U 

0~  3 

4  2_ 

-*  ■ 

0  "*  4 

4  3 

RCL 

0'3 

1  4 

A 

14 

076 

*  C 

0" 

-  c 

e 

O': 

2-  c 

- 

0~9 

01 

090 

00 

0 

081 

93 

= 

OS  2 

4  w 

3  '1 

033 

8. 2 

A  A 

084 

C  P 

083 

04 

«4 

086 

1  w 

*  ^ 

08' 

•#  !• 

h  *  w 

033 

1  *• 

14 

089 

6  7 

w  J 

090 

01 
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COUNTS 


PROGRAM 

STEPS 


30-19 
20  -  42 

43  -  62 

63  -  73 

74  -  33 

34  -  91 

92  -  143 


149  -  213 


Search  routine  *or  finding  correct  voltage-coefficient  group. 

Determining  memory  location  f or  a0  coefficient  if  counter 
>,  7  (counter  ■  storage  14). 

Determining  memory  location  for  a0  coefficient  if  counter 
>.  5  Out  <7. 

If  counter  •  6,  go  to  149  to  determine  T  oecause  need  5th- 
order  polynomial 

Determining  memory  location  for  a0  coefficient  if  counter 
<5. 


counter  •  4,  go  to  149  to  determine  T  oecause  need 
5th-order  polynomial. 


Calculating  ’  using  4th-order  polynomial;  storage  32  contains 
location  of  first.  coefficient;  must  increment  82  to  get 
coefficients  of  •'* .  v-,  v-,  v*,  everything  sxmed  into  storage 
33. 


Calculating  T  using  5th-order  polynomial  storage  32  contains 
location  of  first.  a0}  coefficient,  must  increment  32  to  get 
coefficients  of  v',  v*-,  v^,  v  ,  v*.  everything  sjrnnec  into 
storage  33. 


In  steps  53,  34,  145,  215,  "2na  CP"  is  used  to  clear  the  ’-register. 

’he  memory  locations  are  organized  essentially  the  same  as  in  Program  I. 
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Figure  1.  Flow  Chart  For  Program  II 
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Interpolation 
Cryogenics 
Data  Handling 
computational  'ecnnlques 


A  method  has  been  found  to  interpolate  values  *n jm  a  large  table  of  data 
quickly  and  accurately  using  a  hand-held  programmable  calculator.  Inter- 
polatlon  Is  accomplished  by  evaluation  of  polynomials  which  have  previously 
been  fit  to  the  large  data  table,  ’’he  method  is  perfectly  general  and  In 
particular.  Is  Ideally  suited  to  the  interpolation  of  calibration  tables  In 
the  laboratory  while  an  experiment  is  in  progress. 
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